Astronomy & Astrophysics manuscript no. 1267 


© ESO 2009 


May 16, 2009 





Multi-frequency measurements of the NVSS foreground sources in 

the Cosmic Background Imager fields 

I. Data release 

E. Angelakis\ A. Kraus\ A. C. S. Readhead^, J. A. Zensus\ R. Bustos^ "^, T. R Krichbaum\ A. Witzel^ and 

T. J. Pearson^ 

' Max-Planck-Institut fiir Radioastronomie, Auf dem Hiigel 69, 53121 Bonn, Germany 
- California Institute of Technology, 1200 East California Boulevard, California 91 125, Pasadena, USA 
^ Universidad de Concepcion, Casilla 160-C, Concepcion, Chile 
University of Miami, Department of Physics, 1320 Campo Sano Drive, FL 33146, USA 

Received -, -; accepted -, - 

ABSTRACT 

Context. We present the results of the flux density measurements at 4.85 GHz and 10.45 GHz of a sample of 5 998 NVSS radio sources 
with the Effelsberg 100 m telescope . 

Aims. The initial motivation was the need to identify the NVSS radio sources that could potentially contribute significant contaminat- 
ing flux in the frequency range at which the Cosmic Background Imager experiment operated. 

Methods. An efficient way to achieve this challenging goal has been to compute the high frequency flux density of those sources by 
extrapolating their radio spectrum. This is determined by the three-point spectral index measured on the basis of the NVSS entry at 
1.4 GHz and the measurements at 4.85 GHz and 10.45 GHz carried out with the 100 m Effelsberg telescope. 

Results. These measurements are important since the targeted sample probes the weak part of the flux density distribution, hence the 
decision to make the data available. 

Conclusions. We present the table with flux density measurements of 3 434 sources that showed no confusion allowing reliable 
measurements, their detection rates, their spectral index distribution and an interpretation which explains satisfactorily the observed 
uncertainties. 

Key words. Radio continuum: general - Catalogs - Galaxies: active - cosmic microwave background 

1. Introduction 1.1. The CMB contaminants 



; Targeted multi-frequency surveys can be very efficient in serv- ^^e observer, the CMB is subject to the influence of numer- 
- ing several fields of astrophysical research such as revealing ^us sources of secondary brightness temperature fluctuations, 
new GPS and HFP sources, estimating higher frequency source cumulatively referred to as foregrounds. The reliability of the in- 
counts from extrapolating the radio spectra and hence com- formation extracted from the study of the primordial fluctuation 
^ , puting the confusion limits etc. Consequently, they can be of ^^^^^ spectrum is tightiy bound to how carefully such factors 
• essential importance in the study of the Cosmic Microwave j^^^g been accounted for 
Background radiation (CMB) through the characterization of ^he potential contaminants can crudely be classified in 
the foregrounds. Here, we present the results of the study of a jj^^se of galactic an d tho se of extragala ctic origin (for a re- 
sample of 5 998 NRAO VLA Sky Survey (NVSS, [ Condon et al J ^^^^^ seel^lilllll iTegmark et al.ll2000l) . Moreover, de- 
11998D sources at three frequencies. 1.4 GHz is provided by the p^^^j^g jj^^j^ character, they influence the power spectrum 
NVSS catalog and 4.85 GHz and 10.45 GHz were observed with different angular scales. Galactic foregrounds could be flie 
the Eff^elsberg 100 m radio telescope. The measurements were ^-^^^^ synchrotron emission (for a review, see IsJ 119991) 
initially motivated by the need to estimate the emission that attributed to galactic relativistic electrons, the free-free emis- 
they could contribute at 31 GHz. This is the band in which the ^-^^ originating in Hii regions and flie dust emission due to 
Cosmic Background Imager (CBI, | Padin et alj | 2001|) operates, ^^^^ ^^^^^ ^1^^ interstellar medium that radiate in a black 
as explained below. In a fuhire pubUcation we plan to use the ^^^y ^^^^^^ Extragalactic foregrounds could be the thermal 
extrapolated flux densities in order to compute the source counts and the kinematic Sunyaev-Z el'dovich effect in galaxy clusters 
and the confusion hmits at higher frequencies and compare the jSunvaev & Zeldovichll l970h . manifested through the distortion 
results with those from ottier surveys. l^j^ck body CMB spectrum induced by either hot ionized 

gas in the cluster (in the former case), or matter fluctuations (in 
the latter case). 

In the latter class, radio galaxies and quasars cumulatively 

Send offprint requests to: E. Angelakis referred to as point radio sources, populating the entire radio 



Table 1. 

fields. 



Angelakis et al.: Multi-frequency measurements of the NVSS foreground sources in the CBI fields 
The coordinates of the points defining the targeted 
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1.2. The Cosmic Background Imager 

The CBI is a 13-element planar synthesi s array operating in 10 
1-GHz channels between 26 and 36 GHz jPadin et al.l2001h . It is 
located at an altitude of roughly 5 080 m near Cerro Chajnantor 
in the Atacama desert (northern Chilean Andes). Its task was to 
study the primordial anis otropies at angula r scales from 5'to 
0.5°(400 < i < 3 500, iPadin et alj [20021) . The observations 
of the primordial anisotropies are made in four dis tinct parts of 
the s ky, separated from one another by 6 hr in RA dMason et al.l 
12003; Pearson et al. 2003i table[T]i. 

From the NVSS catalog, it is known that within the CBI 
fields, there are in total 5 998 discrete radio sources with 5 1.4 > 
3.4 mJy. Inevitably, they comprise the potential contaminants 
that may impose second ary fluctuations in the observed back- 
ground temperature field jReadhead et al.ll2b04l) . 



1.3. The Solution 

Instead of removing all potentially contaminated pixels from the 
CMB maps (which would unavoidably cause a significant data 
loss), it would suffice to identify the sources that contribute neg- 
ligible flux density at higher frequencies (below the few-mJy 
threshold) and ignore them during the CMB data analysis. On 
the basis of the assumptions that: 

1. the radio spectrum is described by a simple power law of the 
form S ~ v" (with a hereafter being the spectral index) 

2. the spectrum is not time variable 

this identification can in principle be done by the extrapolation 
of the radio spectrum as obtained at lower frequencies. The ra- 
dio spectrum consists of the flux density at 1 .4 GHz as extracted 
from the NVSS and those at 4.85 and 10.45 GHz as measured 
with the Effelsberg telescope. 



1.4. The Sample 

The list of targeted sources includes all 5 998 NVSS sources 
present in the CBI target fields displaying 5 1.4 > 3.4 mJy. At 
this limit the NVSS is characterized by 99% completeness. The 
sources are distributed in four sky regions - rectangular in ra- 
dec space - with their coordinates shown in table [T] The crite- 
rion for the selection of these regions has been their foreground 



Each point is given by its J2000 (ra, dec) with right ascension (ra) 
measured in hours and the declination (dec) in degrees. 



sky, comprise by definition the most severe contaminant affect- 
ing small angular scales. The sample studied in the current work 
is exactly the NVSS point radio sources that lie within the fields 
targeted by the CBI experiment. 



100 

S,j(mJy) 



Fig. 1. The NVSS 1.4 GHz flux density distribution of our sam- 
ple. Roughly 80% is below 20mJy. This plot demonstrates 
clearly the choice made of radio "quiet" sky regions. 



emission. It has been requ ired that they have IR AS 100/im emis- 
sion less than 1 MJy sr ' (iPearson et al.l2003h . low galactic syn- 
chrotron emission and no point sources brighter than a few hun- 
dred mJy at 1.4 GHz. It is clear therefore that the sample of the 
5 998 sources represents the weak part of the flux density distri- 
bution. This is the most prominent characteristic of the sample. 
In fact, it is readily shown in Fig. [T]that roughly 80% of those 
sources are of 5 1.4 < 20mJy. The large galactic latitudes of the 
targeted fields indicate that the sample is likely to consist com- 
pletely of extragalactic discrete radio sources that is, quasars and 
radio galaxies. 

In the current work we present the results extracted from a 
sample of 3 434 sources which as it is discussed in Sect. [33] and 
14.21 show no confusion from field sources and hence allow reli- 
able measurements. 



2. Observations 

2. 1 . Observing System 

The flux density measurements were conducted with the 
Effelsberg telescope between July 2003 and July 2006. The 
multi-beam heterodyne receivers at 4.85 GHz and 10.45 GHz 
were used. Multi-feed systems use "software beam-switch" 
for removing mostly linear troposphere effects. Each receiver 
was used in two-beam mode (although the 10.45 GHz one is 
equipped with 4 feeds). In both cases, the beams are separated 
in azimuth and each delivers left-handed and right-handed cir- 
cular polarisation channels (LCP and RCP respectively). Both 
systems are mounted in the secondary focus cabin of the 100 m 
telescope. Table[2]gives their characteristics. 

2.2. Observing Technique 

In order to achieve time efficiency, the observations have been 
made with the "on-on" method. Its essence relies on having 
the source in either of the two beams at each observing phase 
whereas the other feed is observing the atmosphere off-source 
(the angular distance of the used feeds is given in table |2|i. The 
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subtraction of the two signals removes linear atmospheric ef- 
fects. For clarity, one complete measurement cycle will hereafter 
be termed as one scan. Each scan, in our case, consists of four 
stages, or sub-scans. 

In order to illustrate the exact observing technique used, we 
label the feeds of any of the two receivers as reference and main. 
Let A be the configuration of having the reference beam on- 
source while the main beam is off-source and B the reciprocal 
case. The telescope is then slewed in such a way as to complete 
a sequence of four sub-scans in a A-B-B-A pattern. Assuming 
then that the system temperature is the same in both feeds for 
any given sub-scan the differentiation of the two signals should 
remove any other contribution than that attributed to the source 
power. The efficiency of the method is demonstrated in Fig. |2l 

Despite its performance , as it is demonstrated in Fig.|2] this 
technique suffers from two major disadvantages: (i) it is sub- 
ject to pointing errors that may result in power loss. This has 
been controlled with frequent pointing checks on strong nearby 
sources. As shown in Sect. 12.31 these errors are negligible; (ii) 
it is subject to cases of confusion i.e. cases of sources that con- 
tribute power to off-source position causing a false subtraction 
result. The solution to that could be either to observe the target 
at a different parallactic angle (at which there would be no con- 
fusing source in the off position), or to correct for it if the power 
of the confusing source is known. This approach is discussed in 
Sect. 1231 

2.3. Logistics 

Thermal noise: For both frequencies, the goal of thermal noise 
(o"rms, see also Sect. I2.4l i around 0.2 mJy (Icr level) has been 
set. Had this been the dominant noise factor, setting a 5cr de- 
tection threshold would allow the detection of sources as weak 
as 1 mJy. The total integration time for achieving this thermal 
noise level is 1 and 4 minutes at 4.85 GHz and 10.45 GHz, re- 
spectively. This time is the cumulative integration time for all 
four sub-scans making up one observing cycle, that is a scan 
(see also Sect. 12.21 ). However, as shown in Sect. 13.21 the domi- 
nant noise factor is the troposphere rather than thermal noise. It 
is shown that the practical limit is of the order of 1 .2 mJy which 
is judged to be adequate. 

Field coverage: A rigid constraint is the minimisation of the 
telescope driving time. This was achieved by driving the tele- 
scope through the field in a "zig-zag" way (travelling salesman 
problem). Each field was organised in stripes parallel to the right 
ascension axis and roughly 0.5 degrees across in declination. 
The sources within such a belt have, in turn, been organised in 
dozens in order of monotonous right ascension change. During 
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Column 1 , is the central frequency while Column 2 is the receiver 
bandwidth. Column 3, gives the system temperature and Column 4 
the sensitivity. Column 5, shows the full width at half maximum and 
Column 6 the angular separation between the two beams. Finally, 
in Colunm 7 we give the polarization channels available. 




Fig. 2. Demonstration of the efficiency of the observing tech- 
nique (upper panel) and a prototype detection profile (lower 
panel) in the case of the 10.45 GHz receiver. Each receiver has 
two feeds, each of which delivers two channels (LCP and RCP), 
giving a total of four channels. Those are shown in the four lower 
panels. The green colour represents the reference horn signal and 
the blue the main horn signal. The left-hand side panels are the 
LCP and the right-hand side panels are the RCP. The plot at the 
top of each panel shows the final profile after subtracting the 
signals from each of the two feeds and averaging over LCP and 
RCP If MR is the RCP of the main horn and ML the LCP in 
the same horn, while RR, RL are for the reference horn, the final 
signal is given by [(ML - RL) + (MR - RR)] /2. It is noteworthy 
that despite the complete absence of even the hint of a source 
in the individual channels (upper panel), after the subtraction a 
clear case of a 22-mK signal (roughly 17 mJy) can be seen. 



an observing session a field would be targeted within hour angle 
range from -3 to 3 hr. This is a compromise between staying in 
one field for as long as possible as well as observing through ac- 
ceptable airmasses (i.e. not very low elevations that would result 
in large opacities). 

Pointing offset minimisation and calibration: For calibration 
purposes, one of the standard calibrators shown in table |3] was 
observed at the beginning and the end of the observation of 
a field, i.e. roughly every six hours. Before the beginning of 
the field, also the focus of the telescope would be optimised. 
Changes in the focal plane within those six hours were accounted 
for by interpolation of the sensitivity factor between the values 
measured at the beginning and the end of the run. To main- 
tain low pointing offsets, cross-scans were frequently performed 
on bright nearby point sources. On average a pointing check 
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Absolute Pointing Offset (arcsec) 



Fig. 3. The distribution of the pointing offsets for the case of the 
10.45 GHz receiver The dashed line represents the offsets in the 
elevation direction while the solid one gives those in azimuth. 
The mean offset is around 3"corresponding to roughly 1 % power 
loss. 



That is, (Tat = |7'ieft-7'iightl/2. Subsequently, the max ((Xrms, ctat) 
is taken as a first estimate of the error in the measurement. In 
Sect. 13.21 we describe how the final errors reported in table H] 
have been calculated. 

2.4.1. Corrections 

Each measurement conducted as described earlier is conse- 
quently subjected to a number of corrections: 

Opacity correction: This process is meant to correct the atten- 
uation of the source signal due to the terrestrial atmosphere. The 
computation of the opacity is done by utilisation of the observed 
system temperatures. 

Elevation dependent gain correction: Small scale diver- 
gences of the primary reflector's geometry from the ideal 
paraboloid lower the sensitivity of the antenna. These gravita- 
tional deformations are a function of elevation with the conse- 
quence of an elevation-dependent antenna gain. The "elevation- 
gain" curve is a second order polynomial of the elevation and 
is constructed experimentally by observing bright point-like 
sources over the entire elevation range. 



was done every 30 minutes to 1.5 hours. This sustained average 
pointing offsets of as low as 3-4"for the 10.45 GHz and 7-8"for 
the 4.85 GHz measurements. These correspond to 4.5% and 5% 
of the FWHM at 10.45 GHz and 4.85 GHz respectively and re- 
sult in a negligible power loss of the order of 1 %. As an example, 
Fig.[3]shows the distribution of pointing offsets for the high fre- 
quency observations. 

2.4. Data Reduction 

Before any further discussion it must be clarified that despite the 
fact that the receivers deliver two circular polarization channels 
(namely LCP and RCP, see Sect. 12.11 ). the possible circular po- 
larization has been neglected with the LCP and RCP channels 
being averaged (see Fig. |2] and appendix iBl). This is a reason- 
able assumption provided that the average degree of circular po- 
larization of these sources is expected to be low (< 0.1%, e.g. 
LWeiler & de Pater 1983; Komesaroff^et al. 1984). 

Figure|2]illustrates the "detection pattern". From that picture 
it is clear that a measurement is the difference between the av- 
erage antenna temperature of the first and the second sub-scans 
(Tieft) as well as that between the third and the fourth (Tright). 
These two differences essentially provide two independent mea- 
surements of the target source. Ideally, the results should be 
identical. Differences should be attributed to atmospheric fluc- 
tuations, given that the overlap of the "off" and the "on" beam 
is not precisely 100%, as well as confusion (field sources con- 
tributing power in the off-beam position). This effect however, 
comprises the most severe uncertainty in the measurement. A 
detailed discussion is given in appendix lAl 

Throughout the data reduction process two types of errors 
are computed. The first, denoted by cTims, is the result of the for- 
mal error propagation (assuming Gaussian statistics) of the data 
scatter around the average (error in mean), is chiefly a property 
of the detector and is practically computed by the radiometer for- 
mula. The second is root mean square (rms) in the antenna tem- 
perature as is measured from the first subtraction (sub-scans 1 
and 2) and that from the second subtraction (sub-scans 3 and 4). 



Sensitivity correction: This process is essentially the transla- 
tion of the antenna temperature to Jy. That is done by observ- 
ing standard calibrators (table [3). Given a source of known flux 
density 5cai [Jy] and measured antenna temperature Ja [K], the 
sensitivity factor F will then be F = Ta/SchI- However, the sen- 
sitivity factors obtained this way depend on the quality of axial 
focusing. This is optimised at the initialisation of a field observa- 
tion. Nevertheless, it can change over the span between two such 
consecutive optimisations (of the order of six hours) and partic- 
ularly when large temperature gradients are present throughout 
the telescope structure. In accounting for that, the sensitivity fac- 
tors have been measured both after the first focus correction (be- 
ginning of the observation) and also before the next focus cor- 
rection (end of the field observation). For an observing instant in 
between, the result of linear interpolation between those two val- 
ues has been used. The flux densit i es of t he calibrators are taken 
from lOtt et"ari(ll994 . lBaars et al.l(ll977l) and Kraus priv. comm. 
It must be noted that apart from NGC 7027 the sources used as 
calibrators are point-like for the beamwidth of Effelsberg tele- 
scope. NGC 7027 on the other hand, is extended at 10.45 GHz. 
At this frequency its size is roughly 9.0xl2.0"(the beamwidth 
at 10.45 GHz is ~ 67"). Nevertheless, the power loss due to this 
effect is still less than 1 % and therefore, no beam correction is 
necessary. 

Confusion: A potential limitation for any observation is con- 
fusion which has been well studied since the early days of 
radio surveys (IScheuej|1957h . Put simply, it refers to blends 
of unresolved sources that build up significant flux densities. 
Traditionally, confusion has been treated statistically in terms of 
expected flux density per unit area for a given frequency and for 
modern observing facilities it often constitutes a factor imposing 
more severe limits than the thermal noise itself. In the case of 
the currently discussed work, the problem becomes even more 
severe because of the beam switch technique used. In this case, 
any combination of field NVSS sources can be in the vicinity of 
the targeted source within the "on" or any of the "off" positions. 
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3. Errors 

In general, the requirement of time efficiency can be in conflict 
with measurement accuracy by limiting, for instance, the time 
invested in calibration. A careful and realistic quantification of 
the involved uncertainties is therefore essential. The following 
discussion deals with the system repeatability study which in 
fact sets the pragmatic limit to the reliably detectable flux den- 
sity. 




(dj 


(or 








(Ol. 



Fig. 4. Confusion examples. The left-hand column shows the de- 
tection profiles whereas the right-hand one shows the NVSS en- 
vironment of each target source. There - assuming a Gaussian 
beam pattern - the solid line marks the 50% beam gain level 
while the dashed one denotes the 10% level. The red circles cor- 
respond to the "on" positions and blue ones to the "off" posi- 
tions. The target sources are shown in red and the environment 
NVSS ones in black. The left-hand side plots show the result of 
the differentiation with respect to the strength of the "confusers" 
and their position relative to the centre of the beam. 



That can severely affect the differentiation algorithm by contam- 
inating the subtracted signal. Some typical confusion cases are 
shown in Fig. |4] The confusion status has been monitored for 
every observed scan on the basis of the NVSS positions and 
has been corrected afterwards whenever possible (see descrip- 
tion appendix IBTi. 



Table 3. The flux densities and spectral indices of the standard 
calibrators. 
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NGC 7027 


5.40 


5.46 


5.83 


5.96 


+ 1.12 
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Columns 2 and 4 give the flux densities assumed for these sources. 
The flux densities and the spectral indices in columns 3, 5 and 6, 7 
respectively, are the ones averaged over our measurements, 
t The average has been done over data with SNR > 5cr. 



3.1. System repeatability 

Given the goal of reaching the telescope's theoretically expected 
least detectable flux density, it is crucial to estimate the repeata- 
bility of a measurement. Let the term "observing system" collec- 
tively describe everything but the target source. Hence, it refers 
to the combination of the telescope, the thermal noise, the at- 
mosphere, the confusion etc. An ideal observing system should 
output exactly the same result for the flux of a source indepen- 
dently of the number of repetitions of the measurements, as long 
as the source itself is not variable. If we therefore assume that 
the source is non-variable, the variance of its measured flux den- 
sity over several repetitions can be perceived as system variabil- 
ity caused by any combination of the possible factors referred 
to previously. The estimation of the mean variance of the sys- 
tem as a whole sets the lower limit in the detectable flux density. 
Considerable observing time has been spent in monitoring ex- 
actly this property of the system. 

A number of sources, hereafter called the "repeaters", have 
been selected to be observed during every observing run. They 
have been chosen to satisfy two conditions: 

1. To be intrinsically non-variable. It is known that sources 
of steep spectrum are unlikely to be intrinsically variable. 
Therefore, a number of sources with spectral index steeper 
than around -0.5 were chosen. 

2. To uniformly cover the whole attempted flux density space. 
This is essential as we expect that the system repeatability, 
as defined above (the rms in the repeatedly measured flux 
density of an intrinsically non-variable source), is a function 
of the flux density of the target. 

Roughly 10 sources per field were selected and repeatedly ob- 
served at the beginning of each observing run of the respective 
field. These sources are included in table|4]along with their aver- 
age fluxes at 1.4, 4.85 GHz and 10.45 GHz as wefl as their low- 
and high-frequency spectral indices. As it is shown there, their 
flux densities cover the range up to a few hundred mJy. In or- 
der to extend the flux density range the pointing sources and the 
main calibrators were also used in the analysis (see tables [3] and 
0. 



3.2. Repeatability plots and error budget 



In Sect. l2.4l it was explained that as a first estimate of the error in 
a measurement, has been taken the maximum between the error 
in the mean after the formal error propagation, o-rms and the part 
influenced by the atmospheric fluctuations and confusion, 0-^7- 
(i.e. mflx(cri-n,s, 0-^7-)). The former is a parameter of the detector 
and is not expected to vary significantly. The latter on the other 
hand can vary even for the same target source as a function of the 
atmospheric conditions and the geometry of the dual-beam sys- 
tem with respect to the target source and its NVSS environment 
(confusion). 
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Table 4. The sources used for pointing correction and the "repeaters". 
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1040 
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2 


0.4 
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-1.19 
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30 


3 


-0.60 


-0.67 


084601-0040 


30 


11 


1 


10 


0.4 


-0.81 


-0.53 


084709-0047 


62 


23 


2 


8 


2 


-0.82 


-1.02 


084721-0025 


72 


19 


2 


8 


1 


-1.09 


-1.10 


084840-0034 


131 


37 


1 


17 


2 


-1.03 


-1.03 


084950-0010 


40 


14 


1 


8 


3 


-0.87 


-0.82 


085255-0023 


32 


13 


1 


7 


1 


-0.74 


-0.75 


085418-0036 


58 


17 


2 


9 


2 


-0.98 


-0.92 


144043+0017 


70 


20 


1 


8 


2 


-0.10 


-1.10 


144119+0025 


84 


24 


2 


11 


2 


-1.02 


-1.01 


144232+0019 


30 


9 


1 


5 


1 


-1.02 


-0.96 


144615+0009 


58 


18 


2 


8 


2 


-0.95 


-0.98 


145004+0024 


51 


13 


1 


5 


1 


-1.10 


-1.15 


145421-0016 


84 


34 


1 


17 


2 


-0.73 


-0.80 


145430-0030 


24 


11 


1 


6 


2 


-0.60 


-0.70 


145548-0037 


67 


36 


2 


19 


2 


-0.49 


-0.63 


145554-0037 


32 


30 


2 


13 


1 


-0.07 


-0.44 


204952-0245 


115 


29 


2 


9 


2 


-1.10 


-1.27 


205001-0249 


261 


95 


2 


46 


2 


-0.82 


-0.87 


205041-0249 


26 


20 


2 


16 


4 


-0.20 


-0.27 


205240-0156 


44 


14 


1 


6 


2 


-0.92 


-1.01 


205546-0204 


94 


24 


1 


10 


2 


-1.08 


-1.14 


205612-0206 


90 


22 


3 


12 


2 


-1.15 


-1.01 


205616-0155 


63 


16 


2 


7 


2 


-1.10 


-1.13 



The second column shows their NVSS flux density. The third and the fifth columns give the average flux densities at 4.85 GHz and 10.45 GHz 
measured at EfFelsberg, respectively. Columns 4 and 6 give the rms scatter in those measurements. 

* The average values are produced by pure detection cases. Only measurements of SNR > 5cr are accepted. 

* Source showing intense variability. Marked with red in the repeatability plots at 10.45 GHz. 



A way to statistically quantify the uncertainty in a measure- 
ment including collectively all possible factors of uncertainty, 
is to investigate how well the measurement of a target source, 
assumed intrinsically non-variable, repeats over several obser- 
vations (see Sect. 13. 11 1. For a given frequency, the measure of 
the system repeatability is the rms in the average flux for every 
repeater as a function of its average flux density, cr,-ep (S). The 
associated plots and are shown in Fig.|5] 



The rms flux density crjep (S), can be written as a function 
of the mean flux density S, being the Pythagorean sum of (i) 
the flux density independent term, ctq and (ii) the flux density 
dependent term, m ■ S An particular, it is described by: 

crrep(5)= yJcrl + im-Sf (1) 

Fitting this function to the 4.85 GHz and 10.45 GHz mea- 
surements, has resulted in the parameters in table |5] From those 
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3.3. Confusion flavors 



1000 10000 



Fig. 5. The repeatability plots. The upper plot corresponds to 
4.85 GHz and the lower one to 10.45 GHz. The parameters of the 
fitted curves are given in table |5] In the lower panel (10.45 GHz) 
the red points correspond to sources that are known to exhibit 
variability characteristics and have been excluded during the fit- 
ting procedure (namely, 025515+0037 and 024104-0815 in or- 
der of flux density). 



fits one can readily estimate the minimum detectable flux den- 
sity at each frequency. Setting the detection threshold at 5 cr, the 
smallest detectable flux density is roughly 6 mJy. In appendix lAl 
we discuss the comparison of the fitted values of ctq with the 
error of an individual measurement and give a quasi-empirical 
interpretation of the measured parameters. From the discussion 
previously and in Sect. 12.41 the most reasonable (and rather con- 
servative) estimate of the errors, would be: 



err = max ((Trms, CTAT, (Trep) 



(2) 



This definition is used to derive the final errors reported in ta- 
ble[8] 



Table 5. The fitted parameters for the repeatability curves. 



Frequency 




Error 


m 


Error 


(GHz) 


(mJy) 


(mJy) 


(%) 


(%) 


4.85 


1.2 


0.2 


1.3 


0.02 


10.45 


1.3 


0.1 


1.6 


0.04 



The cTo parameter is the one determining the least detectable flux. 



Depending on the configuration of the dual-beam system in the 
sky relative to the target source and the instantaneous spatial dis- 
tribution of NVSS field sources, a given scan can show different 
confusion "flavors". On this basis, there can be three scan classes 
discriminated by their confusion status at a given observing fre- 
quency: 

1 . Clean scans. Those are the measurements during which there 
were no contaminating sources within any of the beam posi- 
tions (see top panel in Fig.|6l). For these cases, further action 
need not be taken. 

2. Clustered scans. These are scans on sources that are accom- 
panied by neighboring sources within a radius smaller than 
the associated beam-width (hence the term cluster). These 
sources cannot be discriminated (see middle panel in Fig.|6]) 
and reliable measurement is impossible. For a given fre- 
quency only an instrument of larger aperture could resolve 
them (e.g. interferometer). For this reason, these scans are 
absent from the discussions in the current paper. 

3. Confused scans. This refers to the case of having any com- 
bination of field sources within any beam (see lower panel 
in Fig. |6]l. For these cases one must either conduct a mea- 
surement at different parallactic angle such that the confus- 
ing source will not coincide with an "off" position, or recon- 
struct their flux from the exact knowledge of the flux of the 
"confusers" (see appendixlBt. 

It is important to underline that this effect refers to confusion 
from field NVSS sources alone and not to blends of unresolved 
background radio sources which may contribute significant flux. 

In table |6] we show the detected confusion flavors for each 
field and frequency. From this table it is readily noticeable 
that the confusion becomes less important with increasing fre- 
quency. For instance, the fraction of sources that suffer neither 
from clustering nor from confusion effects increases from 59% 
at 4.85 GHz to 92% at 10.45 GHz. This is easily interpretable 
in terms of smaller beam-width (67"as opposed to 146"at 
4.85 GHz). In fact, considering that the majority of sources show 
steep radio spectra (see Sect. 14. 3t , it is expected that in practice 
significantly fewer sources will suffer from confusion simply be- 
cause their field sources are too weak already at 4.85 GHz and 
10.45 GHz. It is important to state that in the following studies 
we consider only a sub-sample of 3 434 sources which are either 
clean or have been de-confused. 



4. Results 

4.1. Detection Rates 

The essence of our task is identifying the detection rates at each 
observing frequency. Assuming that the detectability of a tar- 
get source is solely due to its spectral behavior, the detection 
rates can reveal the subset of sources that exhibit flat or inverted 
spectra which can adversely affect CBI data (i.e. a > -0.5 with 
S oc v"). The current sub-section deals with this problem. That 
is, essentially counting the sources that have been detected at 
each frequency. For both frequencies the detection threshold has 
been set to 5 cr, with cr being the error in the individual measure- 
ment as defined by Eq.[2] 

A supervisory way to describe the detection rates is using the 
2-bit binary detection descriptor as in table |7] That is, a two-bit 
binary in which the left-hand side bit describes the detection at 
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SOURCE: 024632-0513 50AN: 6965 FREO: 4.65 GH; AZIi 1S25 ELV: 32 7 PA: -10.6 



Table 6. The frequencies of confusion flavors of the observed 
scans (measurements) for each field and observing frequency. 



+ 024829-0506 



[ +024632-0513 



+ 024853-0517 



+ 024BJ4)298Z0-0520 



42.2 42.15 42.1 42.05 42 

RA Ideg) 



SOURCE; 025339-0509 SCAN: 9004 FREO; 4.85 GHz AZI: 169.5 ELV; 33.6 PA: -6,6 



+ 025355-^)5101- 025339-05D9 
V-f 025341-^0510 



+ 025341-0513 



Field 


Clean 


Cluster* 


Confused de-confused 




(%) 


(%) 


(%) 


(%) 






4.85 GHz {FWHM « 145") 




02^ 


59 


18 


23 


4(17) 


08*" 


57 


18 


25 


9(36) 




60 


20 


20 


4(20) 


20'' 


59 


19 


22 


6(27) 


Average 


59 


19 


22 


6(27) 






10.45 GHz {FWHM « 67") 




02'" 


92 


1 


7 


0(0) 


08'" 


93 


2 


5 


0(0) 


14h 


90 


1 


9 


2(22) 


20'" 


92 


2 


6 


1(17) 


Average 


92 


~ 1 


~ 7 


~ 2 (20) 



It must be made clear that this result is based solely on the position 
of the dual-beam system in the sky and the spatial distribution of the 
field NVSS radio sources. In column 5 we give the fraction of the 
total number of scans that have been de-confused and in parenthesis 
the same number but as a fraction of the confused scans which is 
given in column 4. 

Here the term cluster is meant to represent both, the cases of pure 
clustering flavour and those that are clustered and confused simul- 
taneously. 



45 5 45 45 



45 3 45 25 



Table 7. The detection rates. 



SOURCE; 025200-0527 SCAN; 8999 FREQ; 4.85 GHz AZI: 157.2 ELV; 33.1 PA: -7.9 



-10% sensitivity level 



0E5250-0526 



l' ^^^^^^^^4^^ 



+ ra5127-052t 



43 2 43 15 



43.05 45 42 95 42.0 42.85 

RA (deg) 



Field 


sample* 


00 


10 


11 


01 






(%) 


(%) 


(%) 


(%) 


02" 


914 


70.5 


16.0 


12.7 


0.8 


08'" 


692 


66.8 


19.2 


12.4 


1.6 




923 


70.7 


18.2 


10.6 


0.5 


20'' 


905 


67.8 


19.6 


11.3 


1.3 


Total 


3434 


69.0 


18.3 


11.7 


1.0 



The detection threshold for either frequency has been set to 5 cr 
with sigma being the error in the individual measurement as given 
by Eq.|2] 

The sample includes measurements that are clean or de-confused. 
The de-confusion includes also the rare cases of having the source 
observed at different parallactic angle. 



Fig. 6. Confusion flavors. From top to bottom: a clean, a cluster 
and a confusion case. The notation is identical to that in Fig. |4] 



the low frequency and the one on the right-hand side that at the 
high frequency with "0" denoting a non-detection and "1" de- 
noting a detection. From all the sources in the sample we have 
selected only those that are either clean at 4.85 GHz or have been 
de-confused as described in appendix IB] Those sources must 
then also be clean at 10.45 where the beam- width is significantly 
smaller. 

4.2. The measured flux densities 

In table [8] available at the CDS, we summarise the acquired 
Effelsberg measurements along with the computed spectral in- 



dices for each source. For the construction of this table, only 
clean or de-confused cases have been considered. In that table. 
Column 1 lists the name of the source. Columns 2 and 3 give 
the NVSS flux density and its error respectively. Columns 4 and 
5 give the flux density at 4.85 GHz and its error respectively. 
Columns 6 and 7 list the flux density at 10.45 GHz and its error. 
The 4.85 GHz and 10.45 GHz have been measured with 100 m 
radio telescope in Effelsberg. Columns 8 and 9 give the low 
frequency spectral index cKj ^^ between 1.4 GHz and 4.85 GHz 
and its error. Similarly, Columns 10 and 11 give the high fre- 
quency spectral index a^^^'^^ between 4.85 GHz and 10.45 GHz 
and its error Finally, Columns 12 and 13, give the least-square 
fit spectral index a, from the 1.4 GHz (NVSS), 4.85 GHz and 
10.45 GHz data points and its error. 

As mentioned earlier, a measurement is regarded to be a de- 
tection only if it is characterized by SNR > 5. Whenever this 
is not the case an upper limit of 5cr is put, with sigma being 
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the error computed as described in Sect. 12.41 In such cases, the 
associated spectral indices are not quoted in the table. 

According to our discussion in Sect. 13.21 the errors quoted 
in table |8] can not be < 1.2 or < 1.3 at 4.85 or 10.45 GHz, re- 
spectively. Yet, there exist entries that the given error is less than 
those limiting values. The reason for this is that in the rare cases 
that more than one measurements of the same source are avail- 
able a weighted averaged is computed to be the flux density if 
the source. Then the associated error may appear smaller that 
the values of 1.2 and 1.3. 

In that table are included all the measurements that are char- 
acterized as clean or they have been de-confused. Clustered 
sources or cases suffering from confusion are not included. It 
must be noted that concerning the CMB experiments these cases 
still provide useful information. Typically, they are characterized 
by lower angular resolutions and hence clustered sources can be 
treated as individual objects. All in all, the sources included there 
amount to about 57% of the whole sample amounting to 3 434 
entries. 

4.3. Spectral Indices 

The motivation for the current program has been, as discussed 
earlier, the estimation of the extrapolated flux to be expected at 
higher frequency bands performed on the basis of the three-point 
spectral index. Here we summarize the findings of the spectral 
indices study. Hereafter, it is assumed that S oc v". 

To begin with. Fig. |7] shows the spectral index distributions 
for the spectral indices in table [8] In particular, the distributions 
of a'l'l^, o;4°g5^ and least-squares fit three-point a are shown. 
All three of those are constructed only with 5 cr data. For com- 
puting the three-point spectral index, an implementation of the 
nonhnear least-sq uares (NLS) Marquardt-Levenberg algorithm 
dMarquardd 1 9631) was used. That imposes natural weighting (i.e. 
l/cr2). 

The median spectral index a is around -0.71 whereas the 
average value is roughly -0.59 indicating the skewness of the 
distribution. On the other hand, q-^^^ shows a median value of 
roughly -0.69 whereas Q'4*gg^ has also a median of -0.75. Of the 
402 sources detected at both frequencies, 136 (34%) appear with 
a spectral index a > -0.5 which implies that the majority of the 
sources appear to have steep spectrum {a < -0.5). Moreover, 
it nicely explains the large percentage of sources w 87% with 
2-bit binary detection descriptor of 00 or 10. What is important 
about this population is that these are the sources that need not 
be "vetoed out" during the CMB data analysis since they are 
not bright enough to contribute detectable flux at the frequencies 
near 30 GHz at which experiments like CBI operate. 

All the measurements have been conducted in an approxi- 
mately quasi-simultaneous way. The coherence time varies be- 
tween hours to days. It is therefore important to contemplate on 
how the lack of simultaneity influences the results. Provided that 
most of the sources follow a steep spectrum trend and steep spec- 
trum sources are not expected to vary significantly, it is reason- 
able to assume that statistically it will be insignificant. 

5. Conclusions 

1 . The applied observing technique has been chosen to be ef- 
ficient in terms of time. Its combination with the beam- 
switch allows a remarkably efficient removal of linear at- 
mospheric effects. However, it suffers from "analytic" con- 
fusion (caused by sources of positions known from other 



least-squares fit a 




Spectral Index 



Fig. 7. The normalized spectral index distributions. The grey 
area histogram shows the distribution of the three-point least- 
square-fit spectral index, a, for sources that have been detected at 
both frequencies with sigma-to-noise of at least 5 (402 sources, 
table |7]i. The same sub-sample is used for the blue line distribu- 
tion which denotes that of the "high" frequency spectral index, 
'^WT- Finally, the red line shows the distribution of the "low" 
frequency index, a'^^'f, for a number of 1 104 sources detected at 
4.85 GHz (see table |7]i. The mean values of a, a^ ^^ and <3'4''g5^ 
are -0.59, -0.54 and -0.69, respectively. 



surveys) as expected. Nevertheless, the confusion effect de- 
creases fast with frequency (from ~ 22% to ~ 6 % between 
4.85 GHz and 10.45 GHz) thanks to the increase of the tele- 
scope angular resolution. Accounting for the shape of the ra- 
dio spectrum would imply a further decrease in the number 
of sources that can actually cause confusion. 

2. We show that for both the 4.85 GHz and the 10.45 GHz 
observations the dominant factor in the smallest reliably 
(SNR > 5) detectable flux density has been the tropospheric 
turbulence. In Appendix |A] we show that the tropospheric 
factor is of the order of 0.9 mJy and 1 .3 mJy for the 4.85 GHz 
and the 10.45 GHz observations, respectively. On the other 
hand, while the second most important factor for the low 
frequency is the confusion caused by blends of unresolved 
sources (see Section 13.21 and Table O, for the higher fre- 
quency thermal "receiver" noise dominates. The confusion 
in the latter case drops dramatically by an order of magni- 
tude to 0.08 mJy due to the smaller beamwidth and the pre- 
sumed spectral behaviour of radio sources. From this dis- 
cussion it is clear that the major limiting factor has been the 
troposphere itself setting a physical limitation in the least de- 
tectable flux density. That appears to be between 5 x 1 .2 = 6 
and 5 X 1.3 = 6.5 for the 4.85 GHz and the 10.45 GHz re- 
spectively. 

3. The agreement between the interpretation/formulation of the 
errors described in Appendix [Aland the observed ones from 
the study of the "repeaters" is noteworthy. 

4. In Appendix |B] an algorithm for achieving "de-confusion" 
is presented. That is, reconstructing a source antenna tem- 
perature on the basis of some elementary presumptions. The 
algorithm has been successfully used in 6 % of the cases in 
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Table 8. The EfFelsberg measured flux densities along with the NVSS ones and the computed spectral indices (it is assumed that 
S ~ v"). 



Source 


S 1.4 


err 


S 4.85 


err 


S 10.45 


err 


4.i!5 
"l4 


err 


tt^gj^ err a err 




(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 






023958+0137 


4.0 


0.5 


<12.5 




<8.0 










024013+0127 


9.4 


0.6 


<6.0 




<6.5 










024031+0140 


3.6 


0.5 


<18.5 




<9.5 










024057+0129 


7.3 


0.6 


<6.0 




<6.5 










024058+0133 


7.5 


0.5 


<7.5 




<12.0 










024119+0126 


7.8 


0.6 


9.3 


1.2 


<16.0 




+0.14 


0.12 




024124+0154 


49.0 


1.5 


15.6 


2.0 


<14.5 




-0.92 


0.11 




024234+0135 


58.2 


1.8 


35.7 


1.3 


24.9 


1.4 


-0.39 


0.04 


-0.47 0.09 -0.41 0.02 


024402+0140 


4.2 


0.5 


<19.5 




<6.5 










024456+0141 


12.3 


0.6 


<7.5 




<8.0 











Cases of SNR< 5 at either of 4.85 or 10.45 GHz are noted with upper limits of 5 cr. The source name is marked with * or ^ in case de-confusion 
correction has been applied at 4.85 or 10.45 GHz, respectively. The complete table is available in electronic form at the CDS via anonymous 
ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ 



the current study and can be easily generalised in projects 
demanding automation. 
5. In Appendix |C] we present an algorithm that is responsi- 
ble for the "quality check" of every observation (a "scan"). 
Incorporating a number of tests can be used for automatically 
detecting cases of bad quality data and can be generalised to 
be used in a "blind" mode. 
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Appendix A: Expected versus observed 
uncertainties 

Here we investigate how the fitted values of ctq and m compare 
with those semi-empirically known. To begin with, the constant 
term ctq in Eq. [T] can be decomposed in three constituents as 
follows: 



'0 



^conf 



r2 
atm 



(A.l) 



cTf^s the thermal noise, computable from the radiometer 
formula 

cTj-onf the confusion err or, known semi-empirically 

(ICondon et al.ll 19891) 
o"atm variable atmospheric emission error, computable 

from the atmospheric opacity change. 

Separately for each frequency these quantities give: 



Frequency 


^rms 


'^conf 


^ atm 




Fitted ctq 


(GHz) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


4.85 


0.16 


0.80 


0.92 


1.23 


1.2 


10.45 


0.22 


0.08 


1.30 


1.32 


1.3 



being satisfactorily close to the expected values. Similarly, the 
flux dependent part of Eq. [T]can be understood as a multi-factor 
effect. Specifically, it can be written that: 



m' = '«poi + "^cal 



(A.2) 



'Wpoi pointing offset error, easily calculable on the basis 
of a Gaussian-like beam pattern and from the mea- 
sured average pointing offsets 4") 

meal instability of noise diode estimated from Intra-day 
Variability experiments (Kraus priv. comm.) 

'«atm variable atmospheric absorption error, estimated 
from the Water Vapor Radiometer data feoy ,2006,) . 

The expected factor m' as compared to the fitted one m, will be: 
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Frequency 








111 


Fitted m 


(GHz) 


(%) 


(%) 


(%) 


(%) 


(%) 


4.85 


0.21 


1.3 


0.004 


1.32 


1.3 


10.45 


1.01 


1.3 


0.005 


1.64 


1.6 





"Main" horn at 


"Reference" horn at 




sub-scans 2 and 3 


sub-scans 2 and 3 


' \ 






f ^ Q • \ 

,1 S3, ,1 
l' l' 
\\ S4. / 


;' \ ; 

S3. ^, 




\ s ^ ' 






"Main" horn at 


"Reference" horn at 




sub-scans 1 and 4 
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Fig.B.l. The horn arrangement as a function of time during the 
execution of an observation. In blue is the horn that is pointing 
off-source each time and in red the one on-source. Within each 
horn there might be a population of confusing sources the flux 
of which is represented by Si, Si' or Si", with i=l, 2, 3, 4 ... S 
in yellow is the target source. The two blue circles on the left are 
misplaced because the sky rotates within a scan. 



The term determining the detection threshold is clearly the 
term cry in Eq. lA.ll From the above discussion, it is clear that 
for the low frequency observations both the atmospheric and 
the confusion terms are significant. However, For the higher fre- 
quency there appears a decrease in the confusion term (due to 
the smaller beam size) and the dominant remaining factor is the 
atmospheric itself. 



Appendix B: Resolving the confusion 

Here we present a method to partially resolve the confusion 
caused by known field sources. The goal is to reconstruct the 
flux density of a target source whenever possible from know- 
ing the parameters of the ones causing the confusion. These pa- 
rameters are known from the observations described here since 
they also belong to the same source sample and hence are ob- 
served. Note that in the majority of the cases (~ 73%, table |6]l 
the de-confusion is not possible. The reason for that may be ei- 
ther that the confusing sources are not detected or they are con- 
fused themselves. The following analysis is entirely done in the 
(RA,DEC) space. 

Assume a certain orientation, for instance of the 4.85 GHz 
dual beam system, with respect to the target source and a dis- 
tribution of confusing sources as shown in Fig. IB. 11 There, the 
target is represented by the yellow star symbol. In that illustra- 
tion there are three distinguishable populations of sources: The 
"on" population, made of sources that lie within a radius of 1 
beam-width (FWHM) about the target source (5, in Fig. IB. lb 
and each contributes antenna temperature Tj. Hence, cumula- 
tively this group will contribute a brightness temperature; 

ToN^^T, (B.l) 
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The "SUB- 1,4" population, which are the sources S'. located 
within a circle of 1 FWHM of the horn position during sub- 
scan 1 or 4. This is the position of the main horn during the P' 
and 4* sub-scan. This population will contribute a brightness 
temperature: 

i 

The "SUB-2 and 3" population of the sources occupying the 
beam position during sub-scans 2 and 3. This is the position 
of the reference horn during those sub-scans. Their contribution 
will then be: 

7^2-3 = 2 f; (B.3) 

i 

In Eq. IB.ll [B.2| and rB.3i ri, Tj or Tj , is the brightness temper- 
ature contribution of a source at the frequency of interest after 
accounting for the distance from the center of the beam. Hence, 
a source of intrinsic brightness temperature T^rc that lies Xq from 
the center of the beam of the 4.85 GHz system, will contribute: 

r = r„c • e"'""'^* (B.4) 

From the above it is clear that Ton, and will be added 
to the system temperature T^ys altering the result of the differen- 
tiation method. 

Preserving the notation used in Sect. 12. 41 it can be shown that 
the "real" source brightness temperature, Tieai can be recovered 
from observable Tieft and Ti-ight, by: 

T'left = T'real + T'oN ^ - (B.5) 

T'right = T'real + T'oN ^ - (B-6) 

Because the terms Tieai and Ton in practice cannot be resolved 
(cluster cases, angular resolution limitation), it is meaningless to 
refer to them separately. This is why we refer to the cluster cases 
separately throughout the text and why we do not include them 
in table [8] For the sake of the following discussion we refer to 
them cumulatively as Tobs = Ti-^ai + Ton- 

This simple method has some weaknesses: 

1. Clustered "confusers": In the above discussion it is pre- 
sumed that the flux densities of the members of a population 
(e.g. S j), are known from the measurement of which target 
was themselves (all the sources we discuss are from of the 
same sample after all and hence have been targeted). This is 
true only if the distance of a pair of sources of the same pop- 
ulation is > FWHM/2. Hence, the above method has been 
applied only in those cases. 

2. Missing "confusers": The confusing sources are searched 
among the the NVSS ones. Hence, sources that are not de- 
tected by the NVSS survey which may become detectable at 
higher frequencies are neglected. 

3. Upper limits: As seen in table[8] often the upper limits in the 
flux density are significant. However, they are not accounted 
for during the de-confusion algorithm. 

4. No corrections applied: The correction discussed in 
Sect. 12.4.11 are not applied for the confusing sources during 
the resolving algorithm. 

5. Inaccurate positions of beams and non-Gaussian beams: 
In all the above it has been assumed that the positions of 
the beams are precisely known and that there are no point- 
ing offsets. Furthermore, the beam pattern is supposed to be 
described by a simple circular Gaussian. 
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Appendix C: Data reduction "pipeline" 

The data volume acquired during the course of the current 
project has been reduced in a pipeline manner Effort has been 
put into developing software beyond the standard data reduction 
packages used in Effelsberg that could assist the observer to re- 
duce the data as automatically as possible at all stages. Here we 
attempt a rough and very brief description of only some of the 
steps followed. Throughout the pipeline, every system parame- 
ter is monitored and recorded. Some details are omitted in this 
description. 



The front-end: The front-end of the pipeline is the point at 
which "counts" (power) from the telescope are piped in the data 
reduction code. The input consists of four power data channels 
two (LCP and RCP) for each horn along with the signal from 
a noise diode of known temperature, for each one of them, i.e. 
eight channels in total. 



RFI mitigation: Before any operation is applied to the signal. 
Radio Frequency Interference (RFI) mitigation takes place. In 
Fig. lC.ll an example is shown. Here, black represents the signal 
before and red the signal after RFI mitigation. In the top panel all 
four channels of the sky signal are shown in terms of "counts". 
A short-lived spike of extremely intense radiation, characteristic 
of RFI, is clearly seen. A routine iteratively measures the RMS 
in that sub-scan and removes the points above a pre-set thresh- 
old. The resulting signal is shown in red. The same procedure 
is followed for the noise diode signal. Finally, the bottom panel 
shows the final detection pattern free of RFI. 



The signal pre-calibration: After the signals have been 
"cleaned" comes the stage of the comparison of each data point 
with the noise diode signal, both being measured in counts. The 
demand for achieving flux densities as low as theoretically pre- 
dicted for the 100 m telescope imposes the necessity of having a 
noise diode signal that ideally should be constant with an RMS 
of no more than a fraction of a percentile. However, often occur- 
ring cross-talk between different channels or other effects, may 
result in intra-scan instabilities (as shown in Fig. \C.2i that may 
distort the detection pattern. The way around this problem has 
been the idea to normalize ("calibrate") the data to the average, 
over the whole scan, diode signal. The default would be a point- 
by-point calibration that may on the other hand significantly dis- 
tort the detection pattern. 





Fig. C.l. Demonstration of the efficiency of the RFI mitigation 
algorithm. Top panel: the sky signal before (black) and after 
(red) RFI mitigation. Lower panel: the final detection profile free 
of RFI. 
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System temperature measurement: Having the data pre- 
calibrated (meaning in terms of antenna temperature), allows 
system temperature measurements. That in turn, allows measur- 
ing the atmospheric opacity for each particular observing ses- 
sion by using the system temperature of each scan. Later in the 
pipeline, this information is used for correcting for the opacity. 



The corrections: Following the previous stage is that of sub- 
tracting the signal of two feeds and the calculation of the antenna 
temperature. Afterwards, the opacity, gain curve and sensitivity 
corrections are applied as described in Sect. 12.4.11 



The quality check: The conceptual end of the pipeline is the 
quality check subject to which has been every single scan. The 



Fig. C.2. Characteristic cases of intra-scan instabilities of the 
noise diode signal. Each column corresponds to a different scan 
and each row to a different channel. The signal is in terms of 
counts. 



term "quality check" wraps up a number of tests imposed on 
each scan. Some of them are: 

1 . The system temperature of each channel is compared to the 
empirically expected one. Flags are raised at excess of 10, 20 
and 30 %. This test serves as an excellent tracer of weather 
effects, system defects etc. 

2. A second test is the RMS and the peak-to-peak variation of 
the data in each sub-scan for each channel separately as well 
as for the final profile. An increase can be caused by ex- 
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treme non-linear atmospheric effects as well as linear slopes 
present in the data. The latter is most often the result of in- 
creasing atmospheric opacity as the source is tracked at low 
elevations. 

3. In order to trace cases that show a clear linear drift as a 
result of increasing opacity, each scan has been sliced into 
four segments. A straight line has consecutively been fitted 
to each segment. A flag is raised when the slope of a segment 
is above some preset value. 

4. It is examined whether the final measurement profile is in- 
verted and if so whether the absolute source flux density sat- 
isfies the detection threshold being set. It is possible in cases 
of confusion that a source in the off position may result in an 
inverted profile. 

5. It is checked whether sensitivity factor (K/Jy) applied to a 
scan agrees with the empirically expected value with a toler- 
ance of 5 %. 



